Principal component analysis of the reionization history from Planck 2015 data 
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The simple assumption of an instantaneous reionization of the Universe may bias estimates of 
cosmological parameters. In this paper a model-independent principal component method for the 
reionization history is applied to give constraints on the cosmological parameters from recent Planck 
2015 data. We find that the Universe are not completely reionized at redshifts z > 8.5 at 95% 
CL. Both the reionization optical depth and the matter fluctuation amplitude are higher than 
but consistent with those obtained in the standard instantaneous reionization scheme. The high 
estimated value of the matter fluctuation amplitude strengthens the tension between Planck CMB 
observations and some astrophysical data, such as cluster counts and weak lensing. The tension can 
significantly be relieved if the neutrino masses are allowed to vary. Thanks to a high scalar spectral 
index, the low-scale spontaneously broken SUSY inflationary model can fit the data well, which is 
marginally disfavored at 95% CL in the Planck analysis. 


PACS numbers: 98.80.Cq 


1. INTRODUCTION 


Accurate measurements of the reionization history of 
the Universe plays an important role in understanding 
the early stages of structure and star formation. Fortu- 
nately, the evolution of the inter-galactic Lyman-alpha 
opacity measured in the spectra of quasars can provide 
valuable information on the reionization history (1. The 
recent measurements suggest that the reionization of the 
inter-galactic medium was largely complete at redshift 
ze 6 (Qj. The rapid decline in the space density of 
Lyman-alpha emitting galaxies in the region z = 6—8 also 
implies a low-redshift reionization process |3]. However, 
the detailed evolution of the electron ionization fraction 
with redshift is currently unknown due to the difficul- 
ties in probing the reionization state of the Universe at 
higher redshifts through the Lyman-alpha opacity. The 
most promising observation in the near future is detec- 
tion of the 21 cm transition of neutral hydrogen. The 21 
cm emission is a direct probe of the large-scale distribu- 
tion of neutral hydrogen during the epoch of reionization, 
and is, therefore, a complementary probe to the Lyman- 
alpha emission at lower redshifts. According to current 
understanding of reionization sources, high energy pho- 
tons may arise from the first star-forming galaxies and 
quasars. The uniqueness of these objects as the sources 
of photons to reionize the Universe is more challenging 
to establish. 


High precision measurements of the cosmic microwave 
background (CMB) anisotropies provide new insights 
into the reionization history. There are two main ef- 
fects of reionization on the CMB angular power spectra. 
The first effect produces a suppression of the acoustic 
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peaks in the CMB angular power spectra through the 
optical depth integrated over the whole reionization his- 
tory. Earlier reionization leads to the larger suppres- 
sion of the acoustic peaks. The second effect changes the 
shape of the large-scale polarization angular power spec- 
tra. This gives a characteristic bump in the spectra on 
scales larger than the horizon size at reionization. The 
position of the bump is proportional to the square root 
of the redshift at which the reionization occurs, while 
the amplitude is proportional to the optical depth Z! 5). 
Therefore, measurements of the large-scale polarization 
angular power spectra can be used to constrain the reion- 
ization history 


Current constraint on the reionization optical depth 
from the Planck temperature data in combination with 
the low-multipole polarization data is 7 = 0.078 + 0.019 
in the base ACDM model [s], which is smaller than 7 = 
0.089 + 0.014 estimated by the WMAP 9-year data id], 
but consistent with that from WMAP polarization mea- 
surements cleaned for dust emission using 353 GHz po- 
larization maps from the high frequency instrument ; 
These results are based on the assumption of an instan- 
taneous reionization. If such a simple assumption is not 
true, the estimated values of cosmological parameters are 
biased. Since there is currently no convincing model of 
the reionization history, it is important and necessary to 
constrain it in a relatively model-independent way. 


Hu and Holder proposed a principal component anal- 
ysis (PCA) of the reionization history to quantity the 
information contained in the large-scale E-mode polar- 
ization i. The principal components are the eigenfunc- 
tions of the Fisher matrix that describes the dependence 
of the polarization angular power spectra on the reioniza- 
tion history. In practice, no more than the first five prin- 
cipal components are needed to parameterize the reion- 
ization history H. This approach has been applied to 
WMAP data and simulated future data [12H17]. It is 
recently claimed that the Harrison-Zel’dovich primordial 
spectrum is consistent with the WMAP 7-year data in 


such a general reionization scenario is, 19}. In this pa- 
per, we perform the principal component analysis of the 
reionization history from recently released Planck 2015 
data and study the impact of the reionization history 
on the estimates of the cosmological parameters. We 
also investigate the bounds on the neutrino masses and 
primordial tensor perturbations in a general reionization 
scheme. 

The paper is organized as follows. In Section [2] we 
describe the principal component method for parameter- 
izing the reionization history that we will adapt in our 
analysis. The results are presented in Section[B] The last 
section is devoted to our conclusions. 


2. MODEL 


Since there is no convincing physical model of the 
reionization history, the instantaneous reionization his- 
tory is usually assumed when investigating cosmologi- 
cal constraints from CMB measurements. Such a one- 
parameter phenomenological model may bias the esti- 
mates of the optical depth and other cosmological pa- 
rameters. In general, the hydrogen ionization fraction 
Xe(z) is a function of redshift. Binning the hydrogen ion- 
ization fraction in redshift bins is a good idea to parame- 
terize the reionization history. However, the degeneracy 
between the binned ionization fraction leads to weak con- 
straints [20}. The principal component analysis provides 
a solution to this problem m. 

The PCA approach is often adopted to reduce the num- 
ber of variables while containing most of information of 
the system. It converts a set of correlated variables into 
a set of linear uncorrelated variables by an orthogonal 
transformation. The principal components are the eigen- 
vectors of the diagonalized matrix, which are picked out 
according to the corresponding eigenvalues. Intuitively, 
the larger the eigenvalue is, the more information the 
corresponding eigenvector gives. 

Following {11} we parameterize the reionization his- 
tory. Consider a binned ionization fraction x-(z;), i € 
{1,2,..., Nz}, with redshift bins of width Az = 0.25 
spanning Zmin < Z < Zmax. The bin width of Az = 0.25 
is sufficiently small so that the final results we obtain 
are independent of the redshift binning m. Here we 
take Zmin = 6 and Zmax = 30. With the definition 
z1 = Zmin + Az and zy, = Zmax — Az so that N, + 1 = 
(Zmax — 2min)/Az. For z > Zmax, the ionization frac- 
tion is set to the residual fraction of recombination, while 
te = 1.0 for 3 < z < Zmin consistent with astrophysical 
observations P| and ze = 1.08 for z < 3 when the Helium 
reionization is taken into account. The principal compo- 
nents of x-(z;) are the eigenfunctions of the following 
Fisher matrix Fij, 
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FIG. 1: Fractional power spectrum response to a delta func- 
tion perturbation of unit amplitude at 6.25 < zi < 29.75. 
Here we have chosen Az = 0.25. 


which describes the dependence of the polarization spec- 
trum CF* on the ionization fraction ze(z;). To calculate 
the Fisher matrix, we should set a fiducial reionization 
history «®4¢(z;) as well as other cosmological parameters. 
The fiducial model parameters are chosen as the esti- 
mated values of cosmological parameters from the Planck 
2015 data including the temperature and polarization 
power spectra for the base ACDM model [8], listed in 
Table[I] We adopt x®4(z;) = 0.1 in the range of redshifts 
6 < z < 30 to get the same optical depth as the instanta- 
nous reionization. Actually our results are not sensitive 
to the fiducial model. Since the effect of x.(z;) on CP” 
of can be neglected on high £, as shown in Fig. [| we 
can cut off £ safely to €max = 100 in the Fisher matrix. 
The perturbation to the fiducial x4 is a delta function 
correlated with the nearby region [7]. 


TABLE I: Fiducial model consistent with results from the 
2015 Planck temperature and polarization spectra for the 
standard ACDM model. 
Qh Qh h T Ase “T Ns 
0.02225 0.1198 0.6727 0.079 1.882 x 107” 0.9645 


The Fisher matrix Fj; is decomposed as 


Nz 
Fij = (Ne +1)? > Sy(z)o; Su(z;) , (2) 

p=1 
where øo? are the inverse eigenvalues and S\,(z) are the 


eigenfunctions that satisfy the orthogonality and com- 
pleteness relations 
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So these eigenfunctions form a set of orthogonal and com- 
plete bases. An arbitrary reionization history can be de- 
composed into a sum of the eigenmodes, 


telz) = 284(z) + $ mS, (2), (5) 


where m, are the amplitudes of the principal compo- 
nents for a particular reionization history. We order the 
eigenmodes so that the smallest oF, has the lowest index 
number pu, starting at u = 1. Since the effect of each 
eigenmode on CP” becomes smaller as u increases, we 
shall pick out the first few eigenmodes. It is enough to 
choose the first 3— 5 eigenmodes shown in Fig. 2] because 


of the existence of the cosmic variance 
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so that the error in CP” from high p eigenmodes is 
smaller than the cosmic variace. 


S, (2) 


FIG. 2: The five lowest-variance principal components of 
Xe(z) over the redshift range of 6.25 < z < 29.75. The vari- 
ances increase with the index p. 


In our analysis, we parameterize the reionizaton his- 
tory as with the five amplitudes m, for u = 1,...,5. 
We consider a spatially-flat ACDM model described by 
a set of cosmological parameters 


{Wh Ach”, Os, As, ns}, 


where Qh? and Qh? are the physical baryon and cold 
dark matter densities relative to the critical density, O, 
is the ratio of the sound horizon to the angular diameter 
distance at decoupling, A, is the amplitude of the pri- 
mordial power spectrum of scalar perturbations and ns is 
the scalar spectral index. Moreover, we consider two one- 
parameter extensions to the ACDM model, ACDM+ Xm, 
and ACDM-+r, where Xm, is the total mass of neutri- 
nos and r is the tensor-to-scalar ratio. When the tensor 
contribution to the CMB angular power spectra is con- 
sidered, the inflationary consistency relation, ng = —r/8, 
is assumed where n; is the tensor spectral index. For 
comparison with the Planck results, r is defined at the 
pivot scale ką = 0.002 Mpc! (denoted ro.002). 
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We modified the Boltzmann CAMB code to ap- 
propriately incorporate the PCA reionization history and 
used the publicly available CosmoMC package to explore 
the parameter space by means of Monte Carlo Markov 
chains technique RA. 

In our analysis we use recently released Planck like- 
lihood code and data, including the Planck low-l like- 
lihood at multipoles 2 < l < 29 and high-l Plik likeli- 
hood at multipoles 1 > 30 based on pseudo-C; estima- 
tors. The former uses foreground-cleaned LFI 70 GHz 
polarization maps together with the temperature map 
obtained from the Planck 30 to 353 GHz channels by the 
Commander component separation algorithm over 94% 
of the sky. The polarization part of this likelihood is 
denoted as “lowP”. The latter uses 100 GHz, 143 GHz, 
and 217 GHz half-mission cross-power spectra, avoiding 
the Galactic plane as well as the brightest point sources 
and the regions where the CO emission is the strongest. 
“Planck TT+lowP” denotes the combination of the TT 
likelihood at l > 30 and a low-l temperature-polarization 
likelihood. “Planck TT,TE,EKE+lowP” denotes the com- 
bination of the likelihood at | > 30 using TT, TE, and 
EE spectra and the low-l temperature-polarization like- 
lihood. 


3. RESULTS 


Table [L| summarizes the constraints on the reioniza- 
tion history and cosmological parameters from the Planck 
TT+lowP and Planck TT,TE,EE+lowP data. The re- 
sults from the temperature and polarization power spec- 
tra are consistent with those from Planck TT-+lowP, 
but with increased precision. Since temperature-to- 
polarization leakage corrections are not considered in the 
2015 Planck TE and EE power spectra [s], in our analysis 
Planck TT+lowP is usually quoted. Actually the reion- 
ization history is not sensitive to the TE and EE power 
spectra at l > 30. 

The mean posterior values and 68% confidence limits 
of m, are listed in Table I] Compared to WMAP data, 
Planck data place strong constraints on the amplitudes 
of the principal components. The marginalized 2D con- 
tours and posterior probability distributions are plotted 
in Figure[7] As expected, the correlations between {m,,} 
are weak and the probability distributions are Gaussian. 
Moreover, we reconstruct the reionization history from 
Monte Carlo chains using the PCA approach in Figure B] 
The blue curve is the mean value of the hydrogen ioniza- 
tion fraction ze(z) in the redshift range 6.25 < z < 29.75, 
with 68% (dark color) and 95% (light color) confidence 
regions. As we can see the Universe is not completely 
reionized at z > 8.5 at 95% CL. This result is consis- 
tent with the estimated value of the reionization redshift 
Zre = 8.8117 derived from Planck TT+lowP-+lensing in 
the instantaneous reionization model. 

The constraint on the optical depth of the PCA reion- 
ization is T = 0.092+0.017, lager than that of the instan- 


TABLE II: 68% constraints on the PCA reionization and other cosmological parameters from 2015 Planck CMB angular power 
spectra. For the total mass of neutrinos and tensor-to-scalar ratio, the 95% upper limits are given. 


Parameters 


ACDM 


0.02223 + 0.00023 
0.1192+0.0021 
1.04092+0.00047 
0.092+0.017 
ns 0.9689+0.0064 


in(10!9 Ag) 3.115+0.033 


Hg (km s~! Mpc™ t) 67.54+0.96 


0.838 + 0.014 


0.3116 + 0.0132 


+0.0625 
0.0569 7 0.0534 

+0.1193 
—0.0457 0.1091 

+0.1450 
0.0498 0.1440 
0.1503 
0.1500 


+0.1566 
0.0524 0.1570 


—0.010f 


taneous reionization obtained in the 2015 Planck analy- 
sis, but with nearly the same uncertainty [g]. The reason 
is that the ionization fraction at high redshift makes a 
larger contribution to the optical depth than at low red- 
shift. The optical depth between any two redshifts z1 
and 22 is an integration of £e, 


(+2) 


z2 
=0. — dz == 
T(z1, 22) = 0.0691 (1 Yp)%h f EN 


21 


e(z) , 


(7) 
where Yp is the helium abundance. Since the Hubble pa- 
rameter H falls off as (1+z)°/? in the matter domination 
epoch, the same ze can contribute more to the optical 
depth at higher redshift. It is known that the amplitude 
of primordial spectrum of scalar perturbations A, degen- 
erates with optical depth 7 in the form A,e~?7 on small 
scale io, which means that a large 7 leads to a large 
As. Hence a higher value of the matter fluctuation am- 
plitude, og = 0.838 + 0.014, is favored in the PCA reion- 
ization than in the instantaneous reionization, as shown 
in the left panel of Figure This means that such a 
general reionization history strengthens the tension be- 
tween Planck CMB observations and some astrophysical 
data, such as cluster counts and weak lensing. As shown 
in the right panel of Figure |4| this tension can signifi- 
cantly be relieved if the neutrino masses are allowed to 
vary, thanks to the degeneracy between the reionization 
history and the neutrino masses. Larger neutrino masses 
lead to a lower og through the effects of neutrino free 
streaming on structure formation. The constraint on the 
total mass of neutrinos can be relaxed to Um, < 1.03 
eV at 95% CL in the PCA reionization scenario. In Fig- 
ure[b] we also show the correlation of the neutrino masses 
with the Hubble constant. The constraint has a broad 
tail to high masses, which illustrates the acoustic scale 
degeneracy with Ho. 
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FIG. 3: Reconstructed reionization history with the PCA 
method from PlanckTT+lowP. The blue curve is the mean 
value of the hydrogen ionization fraction x-(z) in the redshift 
range 6.25 < z < 29.75, with 68% (dark color) and 95% (light 
color) confidence regions. 


It is claimed in that the Harrison-Zel’dovich 
primordial spectrum (i.e., scale-invariant spectrum) of 
scalar perturbations is consistent with the WMAP 7-year 
data in the PCA reionization scenario. However, our 
analysis indicates that the Harrison-Zel’dovich spectrum 
is disfavored by Planck TT+lowP at 4.9 ø CL, although 
the scalar spectral index is shifted towards slightly higher 
values. To test inflationary models the contribution 
of primordial tensor fluctuations to the CMB angular 
power spectra is taken into account. Figure[6]depicts the 
marginalized 2D contour in the ro.o0o2 — ns plane, with 
68% and 95% CL, for the PCA (blue) and instantaneous 
(red) reionization history, derived from Planck TT+lowP. 
We find that the allowed contour is shifted towards to 
higher values of ns, compared to the instantaneous reion- 
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FIG. 4: Marginalized 2D contours in the og — Qm plane, with 68% and 95% CL, for the ACDM model (left panel) and 
ACDM+ Xm, (right panel) with the PCA (blue) and instantaneous (red) reionization history, derived from Planck TT+lowP. 


ma ACDM+5m,[PCA] 
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FIG. 5: Marginalized 2D contour in the Hp —“mz plane, with 
68% and 95% CL, for the ACDM+ mz model with the PCA 
(blue) and instantaneous (red) reionization history, derived 
from Planck TT+lowP. 


ization. This shift is important for the constraints on in- 
flationary models. For example, it is pointed out in 
that the class of inflationary models with a power-law po- 
tential øP for p < 1 and the low-scale spontaneously bro- 
ken SUSY inflationary model are marginally disfavored 
by Planck data at 95% CL. In the PCA reionization sce- 
nario, however, these models are in agreement with re- 
cent CMB observations. Moreover, the constraints on 
the tensor-to-scalar ratio is tighter compared to the in- 
stantaneous reionization because the large amplitude of 
scalar spectrum suppresses the tensor-to-scalar ratio. 
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FIG. 6: Marginalized 2D contour in the ro.o02—ns plane, with 
68% and 95% CL, for the ACDM+r model with the PCA 
(blue) and instantaneous (red) reionization history, derived 
from Planck TT+lowP. 


4. CONCLUSIONS 


We have considered the principal component 
parametrization of the reionization history and studied 
the effects of the PCA reionization on the estimates of 
cosmological parameters using recently released Planck 
data. From the reconstructed reionization history we 
found that the Universe is not completely reionized 
at z > 8.5 at 95% CL. This conclusion can be tested 
by future astrophysical experiments such as SKA. The 
reionization optical depth is higher than but consistent 
with that obtained in the instantaneous reionization 
scheme. The higher value of 7 leads to the higher 
value of As, hence the higher og, which strengthens the 
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FIG. 7: Marginalized 2D contours (68% and 95% CL) and posterior distributions for the amplitudes of the principal components 


my, derived from Planck TT+lowP. 


tension between Planck CMB observations and some 
astrophysical data, such as cluster counts and weak 
lensing. But the tension can significantly be relieved if 
the neutrino masses are allowed to vary. The constraint 
on the total mass of neutrinos from CMB data is relaxed 
due to the degeneracy with the reionization history. 
Thanks to the shift of the scalar spectral index to higher 
values in the PCA approach, the low-scale spontaneously 
broken SUSY inflationary model, which is marginally 
disfavored at 95% confidence level in the Planck analysis, 
can fit the data well. 
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